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Abstract 

In this paper we provide a proof of the multiplicative kinematic description of crystal elasto- 
plasticity in the setting of large deformations, i.e. F = F e F p , for a two dimensional single 
crystal. The proof starts by considering a general configuration at the mesoscopic scale, where 
the dislocations are discrete line defects (points in the two-dimensional description used here) 
and the displacement field can be considered continuous everywhere in the domain except at 
the slip surfaces, over which there is a displacement jump. At such scale, as previously shown 
by two of the authors, there exists unique physically-based definitions of the total deformation 
tensor F and the elastic and plastic tensors F e and F p that do not require the consideration of 
any non-realizable intermediate configuration and do not assume any a priori relation between 
them of the form F = F e F p . This mesoscopic description is then passed to the continuum 
limit via homogenization i.e., by increasing the number of slip surfaces to infinity and reducing 
the lattice parameter to zero. We show for two-dimensional deformations of initially perfect 
single crystals that the classical continuum formulation is recovered in the limit with F = F e F p , 
det F p = 1 and G = Curl F p the dislocation density tensor. 


1 Introduction 


Standard continuum models in the setting of large deformations are based on the kinematic as¬ 
sumption F = F e F p , which decomposes the total deformation gradient F multiplicatively into the 
effective deformation induced by elastic and plastic mechanisms, F e and F p respectively. This 


approach was first introduced i 
many engineering applications ( 

n the 1960’s (Lee and Liu, 196' 

7) and has been very successful in 

Simo, 1988; 

Ortiz and Repetto, 

1999 Nemat-Nasser 

2004; Barton 

et al. 2013; 

Dassault Systemes| 

, 2012). However, the lack of a well-grounded justification for this 


decomposition has led to numerous debates in the literature. The heuristics behind the expression 
F = F e F p is based on the chain rule for the derivative of the deformation mapping F = Dcp, when 
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the elastic deformation of the system can be fully relaxed, leading to a purely plastic distortion. At 
the core of the problem lies the fact that dislocations, when present, physically couple the elastic 
and plastic field around them, preventing the existence of a physical configuration with associated 
deformation F e (or F p ) and thus exclusively elastic (or plastic). This inherent complexity has raised 
many issues, such as the physical meaning of the individual tensors F e or F p , the existence (Casey 
1992; |Deseri and Owen , 2002) and uniqueness (Nemat-Nasser 


and Naghdi 


1981 


Zbib 


Dafalias 


1993 


1987 


Naghdi 


Rice 


1971 


1990 


Green and Naghdi] |1971[ |Casey and Naghdi , 1980 ; Mandel , 1973 


1979 


Lubarda and Lee 


Xiao et al. 


2006) of the decomposition F = F e F p , or the appropriate 


measure of dislocation content in the body (Bilby, 1955; Eshelby. 1956; Kroner, 1960; Fox, 1966 


Willis, 1967; Acharya and Bassani, 2000; Cermelli and Gurtin, 2001) based in the incompatibilities 


of elastic or plastic component of the deformation. Other decompositions have also been proposed 


in the literature, some additive, typically in rate form (Nemat-Nasser 


1994), and others multiplicative (F = F p F e , Clifton (1972), 


Lubarda 


1979 Zbib, 1993 Pantelides 


(1999); F based on the prod¬ 


uct of three tensors, one elastic and two inelastic, Lion (|2000|); Gerken and Dawson (2008); Henann 
and Anand (2009); Clayton et al. ( 2014| )) 


A micromechanical definition of the different tensors F, F e and F p has been recently provided 


(Reina and Conti, 2014) without any a priori assumption on the relationship between them. These 
definitions are exclusively based on kinematic arguments and physical considerations, are uniquely 
determined from the microscopic displacement field, and lead to the expected multiplicative de¬ 
composition in the continuum limit in the absence of dislocations. From the previous definition of 
F p , Curl F p arises as the natural measure for the dislocation density tensor in the material. 

Based on the above micromechanical definitions of F, F e and F p , we provide in this paper a 
rigorous proof of the continuum kinematic assumption F = F e F p , with detF p = 1, for a general 
dislocation and slip structure. Such generality is achieved in this work by considering microscopic 
elastoplastic displacement fields of finite energy that belong to a functional space in which dis¬ 
placements are continuous everywhere in the domain except potentially at the slip surfaces over 
which there is a displacement jump. The continuum limit is then mathematically studied by letting 
the lattice parameter tend to zero and increasing the number of dislocations and slip surfaces to 
infinity. We show, for a two-dimensional deformation of an initially perfect single crystal, that 
the displacement field is indeed continuous in the limit, and that the limit of F coincides with the 
product of the limit of F e and F p . The later statement is highly non-trivial since limits and prod¬ 
ucts do not generally commute. This complexity, combined with the non-linearity induced by the 
finite kinematic description, has so far limited most of the mathematical results of coarse-graining 


of dislocation ensembles to the linearized kinematic setting, see e.g. Conti and Ortiz (2005); Gar- 


roni et al. (2010), where the principle of superposition applies and the decomposition can be easily 
proven to be additive, i.e. e = s e + s p , or to the case of well-separated dislocations (Scardia and 


Zeppieri, 2012 Muller et al., 2015). Even at the continuum scale, few mathematical results exist 


on finite elastoplasticity. Some examples can be found in Mielke and Muller (2006); Mainik and 


Mielke (2009) and Conti et al. (2011). A discrete model which does not require the existence of a 


reference configuration was developed by Luckhaus and Mugnai| Q2010 ). 

The paper is organized as follows. We begin in Section [2] by introducing the mathematical 
notation needed for the remainder of the text and reviewing the microscopic definitions of F, F e 
and F p following the lines of Reina and Conti (2014). Next, in Section [ 5 ] we discuss the passage 
of the aforementioned quantities to the continuum limit from a physical perspective, and provide 
some further guidance on notation in Section [4j This is followed by Sections [5] and [6] where the 
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Figure 1: Elastoplastic deformation (right images) of a perfect crystal (left images) at the atomistic 
(topQ and mesoscopic scale (bottom) induced by the glide of an edge dislocation from the right 
surface of the domain till its center. The surface of displacement discontinuity is represented in red 
in the undeformed configuration and terminates at the dislocation point and the boundary of the 
domain. 



precise setting of the problem and the mathematical results are shown. The paper is then finalized 
with some conclusions in Section [3 


2 Mesoscopic description of elastoplastic deformations 


The present study is exclusively concerned with two-dimensional elastoplastic deformations induced 
by dislocation glide (of edge type) and disregards other processes such as phase transformations, 
diffusion, void nucleation or appearance of microcracks. Mathematically, we describe these de¬ 
formations with a semi-continuous Lagrangian formulation, where the final configuration </?(X) of 
every material point X in the reference domain ft C M 2 , considered here to be a perfect crystal, 
is described via a deformation mapping (p : Q -A cp(fi) C M 2 . In particular, cp is continuous ev¬ 
erywhere in the domain, except at the area swept by the dislocations during their motion (lines in 
two dimensions), over which there is a displacement jump or discontinuity, c.f. Fig. [lj These type 
of deformations belong to the space of special functions of bounded variation or SBV , c.f. Evans 


and Gariepy (1991); Ambrosio et al. (2000); and we further require that cp is one-to-one in order to 


avoid interpenetration of matter. The deformation gradient F is then defined as the distributional 
gradient of <£>, which for SBV functions is a measure of the form 


F = Dtp = \7(p C 2 + |[v?] <g> N n 1 L J. 


( 2 . 1 ) 


The measure F consists of an absolutely continuous part, V<£ C 2 , and a singular part, [[<£>] 0 
N V, ] [j, that has its support over the slip lines or jump set J , where N(X) is the unit normal 

2 Reprinted from Publication ‘Kinematic description of crystal plasticity in the finite kinematic framework: A 
micromechanical understanding of F=FeFp’, Vol 67, Authors C. Reina and S. Conti, Page No. 43, Copyright (2014), 
with permission from Elsevier 
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Figure 2: Atomic positions of a perfect crystal (left image) after a plastic deformation (middle 
image) and an elastoplastic deformation (right image)^] 


at X G J and [<£>] is the jump over the slip line, see below. We denote by C? the Lebesgue 
measure (which measures area in the plane) and by H 1 the one-dimensional Hausdorff measure 
(which measures length). V<£ is called the approximate differential of cp and corresponds to the 
total deformation F at the points where no slip occurs and the deformation is thus purely elastic. 
As a result, \7p> can be physically identified with the elastic deformation tensor 


F e = \7cp. 


( 2 . 2 ) 


The symbol V corresponds to the ‘standard’ gradient where the function is differentiable, in par¬ 
ticular, only outside of the jump set. By definition, F has a vanishing Curl. However, Curl \7p> 
does not vanish in general, and the same occurs for the singular part of the deformation gradient. 


We note that an additive decomposition of the deformation gradient of the form of (2.1), with the 


first term being identified as the elastic deformation tensor, has been previously considered in the 
mechanics literature based on local volumetric averages and application of the divergence theorem 
over continuous regions of the domain (Davison, 1995). 

For its part, [y?](X) = </? + (X) — </? - (X) is the displacement jump at point X in the slip line, 
where the + side is the one indicated by the normal N to the line at point X. Without further 
restrictions on [</?] and J, the class of deformation mappings considered thus far can also model 
processes such as cavitation, crack opening or interpenetration of matter. A careful analysis of the 
kinematics of slip in crystalline materials, c.f. Reina and Conti ([2014), indicates that jump sets in 
two dimensions necessarily consist of an ensemble of straight segments in the reference configuration 
that terminate at dislocation points or exit the domain, c.f. Fig. [TJ and that the displacement jump 
for single slip satisfies, c.f. Fig. [2j 


cp-{X) = cp + {X- b), VX,X-beJ, 
<^+(X) = ¥ >-(x + b), VX,X + beJ, 


(2.3) 


where b is called the Burgers vector and represents the interatomic distance in the direction of slip. 


We remark that the two conditions in (2.3) are equivalent when X, X — b, X + b G J. 


The jump set J and the previously introduced vector b for a single slip system, clearly char¬ 
acterize the plastic deformation with quantities independent of the elastic distortion. Indeed, if we 


4 Reprinted from Publication ‘Kinematic description of crystal plasticity in the finite kinematic framework: A 
micromechanical understanding of F=FeFp’, Vol 67, Authors C. Reina and S. Conti, Page No. 45, Copyright (2014), 
with permission from Elsevier. 
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Figure 3: Sequence of plastic deformations involving two orthogonal slip systems. The zoom in the 
reference configuration is used to label the segments composing the jump set. 


consider for instance a plastic simple shear deformation, such as ip 1 in Fig. [31 the total deformation 
is purely plastic and thus equal to F p , and can be expressed, c.f. Eq. 0> , as 

F = Dp 1 = F p = I C 2 + |bi|ei 0 e 2 H 1 [j. (2.4) 


Similarly, for the composition of two simple shears where slip lines intersect, c.f. Fig.[3]an d[R eina| 
and Conti (2014), the plastic deformation tensor is uniquely given from Eq. ( |2.1| ) as 


F = D (cp 2 ° = F P = 1 £ 2 + |bi| ei®e 2 T^Lji + lb^ e 2 ®ei U l [j 2 + |b 2 |e 2 (g)e 2 (2.5) 


The last term in Eq. (2.5) relates to the kink that results from the pullback of the second slip lines 


into the reference configuration and provides the non-commutative term in the finite kinematic 
formulation. In other words, the pullback is performed in opposite order in which the slip occurred 
and therefore the resulting kink term is dependent on whether the slip is first performed in the 
horizontal or vertical direction. 

The above results can be easily generalized for an arbitrary sequence of N s dislocation-free 
plastic slips, not necessarily homogeneous or belonging to orthogonal slip systems. In that case, <p 
can be written as (p = cp Ng o...o^o...o^ 1 , where (jump set associated to (p u ) is an ensemble 
of parallel straight infinite lines with normal and 

Dp v = I £ 2 + bj, ® N„ V} \_j Vv , h u ■ N„ = 0, b„eR 2 constant on each line. (2.6) 
The total plastic deformation tensor associated to cp can then be written as 


F p = I £ 2 + y] bj 0 Nj- n 1 [jj, (2.7) 

3 

where the total jump set J in the reference configuration consists on the union of straight segments 
Jj (appropriately including all kinks), each of which has an associated Burgers vector hj. In 
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general, as was the case for the composition of two simple shears, b j • N j is not necessarily zero. 
Equation (2.7) remains of course valid if an elastic deformation is superposed to the previous plastic 

o<pP, F p = Dcp p . In such case, by construction, F = Dcp = D (cp e o cp p ) = 

o cp p , and, if F e is sufficiently 


distortion, i.e. cp = cp' 

\7cp e o 

smooth one obtains that for compatible domains F = F e F p + (9(|b| 2 


' e - cp p C 2 + (( p e (cp p +) - cp e (cp p ~)) ® N n 1 [j, F e = V(p = V(p e O <nP ^ ;f l?e 


c.f. Reina and Conti (2014). 


In a general elastoplastic deformation, dislocations are present in the body, coupling the elastic 
and plastic field. They render impossible a global decomposition of the total deformation mapping 
into a purely plastic distortion and a subsequent elastic deformation, i.e. cp ^ cp e o<p p , thus compli¬ 
cating, in principle, the physical definition of the global plastic deformation tensor F p . However, at 
this mesoscopic scale where the dislocations are individually resolved, subdomains away from the 
dislocations are, by construction, defect free and it is therefore possible to obtain in each of them 
a decomposition of the form cp = cp e o cp p and an associated F p given by Dcp p , as in Eq. (2.7). It 
will be shown in Section [5] that F p is uniquely defined from cp regardless of the potentially many 
decompositions of the form cp e o cp p in each of these subdomains, thus providing a global definition 
of F p almost everywhere in the domain (except at the dislocation cores and close to the boundary). 

It can therefore be concluded that F, F e and F p are uniquely defined from the mesoscopic 
deformation mapping cp in general elastoplastic deformations induced by dislocation glide, and are 
of the form 


F = D(p = \7ip £ 2 + <g> Nj n l y., 

j 

F e = Vc p, (2.8) 

FP = I £ 2 + ^b i ®N i U\ Jy 


As will be shown in the following sections, the precise definition of F p inside the dislocation cores 
will be irrelevant in the continuous limit. For concreteness, we keep the definition above with b j at 
each slip line inside the cores as the Burgers vector of the corresponding segment outside of the core 
(this will be made precise in Section [5]) . This choice of F p associates to each segment composing 
the jump set a constant Burgers vector, and as derived analytically in Reina and Conti (2014), 
Curl F p = bk 5x fc , where {X&} is the ensemble of dislocation points and b& the corresponding 
Burgers vector (sum of the Burgers vector of the slip lines terminating at the dislocation point, 
with the appropriate sign). In other words, G = Curl F p exactly measures the dislocation content 
in the body and is therefore an appropriate definition of the dislocation density tensor. 


3 Scaling towards the continuum limit 


Ultimately, we are interested in the description of the elastoplastic deformation in the continuum 
limit. To that regard, we consider sequences of deformations, c.f. Fig. [4j in which each element of 
the sequence is characterized by the parameter e that measures the normalized lattice parameter 
(ratio between the shortest atomic distance and the characteristic length of the body). In particular, 
each deformation mapping will be denoted as (p e , and its jump set and the Burgers vector of each 
segment as J — JT Jj and brespectively. Similarly, the different deformation tensors, as given 
by Eqs. (2.8), will be denoted as F e , F^ and F p . Furthermore, by crystallographic arguments, the 
number of slip systems is finite and bounded by some fixed N s ; I be? | < ne for some fixed n > 0; 
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7V e = 1, e = (7 N e = 3,e = C/3 


N e oo, e = C/A7 0 


Figure 4: Sequence of elastoplastic deformations as the lattice parameter e tends to zero, and the 
number of dislocations N e and length of the jump set \J\ tends to infinity. 



Figure 5: Jump sets corresponding to the same slip system are separated by a distance of at least 

e. In the case of two orthogonal slip systems, the maximum length of slip lines contained in a ball 

2 

of radius r is then proportional to 1 —. 


and planes of the same slip system are separated by a distance that is equal to or larger than e. 
From the latter, it is implied that, c.f. Fig. [5| 

r 2 

n l {Jf\B r (X))<A-~, (3.1) 

for some A that depends on N s . We note that the subindex e is omitted for the jump set J since 
the limiting deformation mapping will be continuous, and therefore J is unequivocally associated 
to an element of the sequence. 

The continuum limit is then attained as e vanishes, and the macroscopic kinematic quantities 
<p, F, F e and F p are uniquely defined as the limit of the corresponding mesoscopic quantities, i.e. 

¥> c ->¥>, F e —» F, F®->F e , F p F p (3.2) 


as e 0 in the appropriate topology. This perspective thus delivers a well defined F e and F p , 
that do not invoque any macroscopic unloading path, which in general does not exist. Only at 
the mesoscopic scale, as discussed in the previous section, a relaxation process (or decomposition 
of the deformation mapping into a physically realizable elastic and plastic deformation) is defined 
everywhere except at the dislocation cores. 

In analogy to a zoom out process in a real material, the number of dislocations N e and the total 
length of the jump set \J\ tend, potentially, to infinity as e —» 0. Equation (3.1) directly implies 

U l {J)e<C\ (3.3) 
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for some (7*, and thus Fj? is at most a finite deformation, as desired. We further assume that the 
number of dislocations N e satisfies 

N e e < C, (3.4) 


for some C. This scaling allows the formation of dislocation walls and prevents the clustering 
of infinitely many dislocations at one macroscopic material point, as is commonly observed in 
experiments. To further support the scaling given by Eq. (3.4), consider an ensemble of identical 
dislocations positioned over a regular square pattern, with distance S between dislocations. A 
simple energetic analysis in the linearized kinematic framework indicates that the system will have 
finite elastic energy only if 5 /yT is bounded, and therefore the total amount of dislocations N e is 
proportional to ^ 

In the following sections, we will show that the macroscopic quantities F, F e and F p defined in 
(3.2) are related in the continuum limit via the standard multiplicative decomposition F = F e F p 
with detF p = 1, that p is continuous in the limit and that G = Curl F p represents the continuum 
dislocation density tensor when expressed in the reference configuration. This is done for general 
sequences of elastoplastic deformations with bounded energy. 


4 Notation 


The notation used in the foregoing will follow, when possible, the standard notation used in con¬ 
tinuum mechanics. In particular, vectors and tensors are denoted with boldface symbols (the order 
of the tensor will be clear from the context), and scalars are written in standard font. Similarly 
dX = dC? denotes the differential area in the reference configuration and dX a differential vector, 
so that x ^ 

[ 2 f(X)-dX= [ f(Xi+t(X 2 -Xi))-(X 2 -Xi)dt. (4.1) 

Jx i J o 

The scalar product between two vectors is denoted as f • g, the full contraction between two tensors 
as A : B, the matrix vector multiplication with no symbol Ag and ® is used to denote the dyadic 
product between two vectors f 0 g. Furthermore, | • | is used to denote the Euclidean norm of the 
scalar, vector or tensor that the symbol encloses. For sets, we use the notation \A\ = T~L l {A ) or 
|A| = £ 2 (A), where the dimension of the set A is clear from the context. 

The spaces that will be used in the following are the space of measures A4, the L p spaces, the 
space of functions of bounded variation BV and the space of special functions of bounded variations 


SBV . Further details on these spaces may be found in Ambrosio et al. (2000), Evans and Gariepy 


(1991) and Adams and Fournier (2003). 


The deformation mappings cp e and cp^ will be considered uniquely defined at the jump set, by 
arbitrarily choosing the value at either side of the jump. With these considerations in mind, both 
mappings are bijective, see Fig. [6j In the compatible regions of the domain, where cp e is of the form 
— the choice of p e , p\ and the domain of <p\ need to be consistent with each other. 

Similar considerations apply for F|(X) = F^ o p v e in the compatible region, where the tilde is used 
to denote quantities in the image of (pe, i.e. X = <^e(X). 

Finally, for easiness in the notation, the subindex k will be reserved for dislocations, j for 
segments of the jump set, and l for domains in the compatible regions away from the dislocation 


cores. 













0e(X • N) 


X-N 


^«X-0 e (X-N)N i 




Figure 6: Example of a single slip plastic deformation (X) = X — </> e (X • N)N^. The value of 
at the jump point is arbitrary, but uniquely defined, so that (fe is bijective. 


5 Setting of the problem 

Let ft C R 2 be an open bounded set with Lipschitz boundary, describing the reference configuration 
of a defect-free crystalline material. For every e > 0, we consider deformation mappings cp e : ft 
<p e (ft) C M 2 , (p € G X e defined below, with finite energy, described by the following functional 


E Te (cp e ) = E e ((p e ) - [ t-ipedX- [ t -y e dU\ 
Jn Jan 


(5.1) 


where f are body forces, t are the prescribed tractions at the boundary (body forces and tractions 
are assumed to be a system of zero net force and momentum for static equilibrium) and 


Ee(<Pe)= [ iVe W e (F e e )dX + a L 

Jn\\JB ce (x k ) J l 


' N e 

U Be 


:(Xfc) 


W e (F e e ) dX + /3\DF e e \(Q), 


(5.2) 


where a, f3 > 0 are material parameters and F® is the elastic deformation tensor. We note that 
only the energy terms associated to E e are needed for the proofs, and that additional terms may 
be added to Eq. 0 without affecting the results. In particular, a term dependent on the plastic 
deformation tensor Fg is physically justified to provide an energetic compromise between elastic 
and plastic mechanisms for achieving a deformation compliant with the external loads. 

The elastic energy functional given by Eq. (5.2) is composed of three terms. The first two 
represent the elastic energy outside and inside the dislocation cores respectively, where each core 
is considered to be a ball of radius ce, around the dislocation points X&. As usual, we consider 
continuous elastic energy densities W e ( Ff) with quadratic growth, i.e., 


1 


— \F e \ z -C L < W e { F e ) < Cu\F e \ 2 + C v , 
El 


(5.3) 


with C L ,Cu > 0. Continuum mechanics is expected to fail in the core, and the parameter a 
provides the necessary rescaling of the dislocation core energy to match its physical value, as given, 
for instance, by atomistic simulations. However, the value of a will be irrelevant for the proofs and 


may be set to one if desired. For its part, the last term of Eq. (5.2) represents the total variation of 


the measure DF e e over the set ft. It is mainly introduced for mathematical reasons, and it is needed 
to obtain F e « F|Fe in the domain, see Proposition 6.1 It can be physically interpreted as a non¬ 


local elastic deformation term that prevents, for instance, the formation of infinitely thin laminated 
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microstructure with alternate elastic deformations. This type of energy terms are common in the 


2009) 


so-called ‘gradient elasticity’ theories (Aifantis, 

The functional space X e , defined below, includes the kinematic restrictions of slip described 
in Section [2j the crystallographic conditions and scalings discussed in Section [3j and considers 
two additional constraints that attend to mathematical simplicity. In particular, dislocations are 
considered to be separated from each other and from the boundary of the domain by a distance of at 
least several atomic spacings, and the sequence of potentially activated plastic slip systems will be 
considered fixed with all slip systems distinct. Although this last assumption may appear artificial 
and indeed does not account for the full set of compatible elastoplastic deformations that a material 
may suffer (a weaving structure common in textiles would be out of this scope), we will show later 
in Lemma |5.6| that an arbitrary plastic deformation tensor may be written as the composition of 
three simple shears with prescribed slip direction and order of activation. The domain is then 
separated into the dislocation cores £> ce (X&), where the Burgers vector is spread over a ball of 
radius me < ce to remove the elastic singularity, and ft \ U^ e L> me (X&), where the deformation 
mapping is compatible. The regions of overlap, £> ce (X&) \ L? me (X&), are locally compatible and 
allow a natural extension of the definition of F p everywhere inside the core region. As will be 
shown in the following, the specific details of the deformation inside each of these cores as well as 
the definition of F p in these regions will be irrelevant in the continuum limit, justifying the use of 
specific functional forms for the deformation in the cores. 

In the following we give a precise definition of the space X e . It depends on three global constants 
independent of e, namely, the number of distinct slip systems N s , the normalized maximum length 
of the Burgers vector for an individual slip |b € /e| < n', and the scaling associated to the number 
of dislocations (7, c.f. Eq. ( |3.4| ). We define 

X e = {(p e G SBV(ft] M 2 ) : <^ € , F®,F P satisfy conditions (C1)-(C4)} (5.4) 

for some global constants n',C,N s > 0 and some N s global unit vectors Ni, ... , N^, ...,Njv s , all 
distinct. The constants c, m, L, n > 0 will be defined below from n', (7, N s . 


(Cl) There are N e dislocation points (Xfc)fc, with 

e 

which describe the dislocations. They are well-separated in the sense that 


(5.5) 


l-X-fci -X-&21 ^ 3ce, Mk\ /t2 5 
dist(Xfc, dft) > 2ce, Mk — 1 ,.. 


.,N e , 


(5.6) 


where c — m-f 2L 2 , with m = 2N s n' and L > 0 as defined later in Lemma 5.1 The value of 


m is discussed in (C3) and the value of c results from Lemma 5.3 


The jump set in each dislocation core X&, J flB ce (Xfc), consists of N^k < X s distinct straight 
segments joining X& with a point in dB ce (Xk ). 

(C2) For every convex domain wCl 2 \ (j£i ^me(Xfe), (f e is a compatible deformation in uo D 14, 
in the sense that there are two maps (f>\ and (p^ such that 


(p € (X) = (cp e e o pf) (X) VX G cj n ft, 


(5.7) 
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where fi\ : R 2 —> R 2 is Lipschitz continuous and one-to-one and : R 2 -G R 2 can be written 
as the composition of N s single slip deformations, fie = fi^ n s ° ••• ° (fie^ ° ... ° where each 
fi^v is of the form 

^(X) = X-^(X- NON^. (5.8) 


Here is defined after (5.4) and each function 0 e?z , : R —>► R is piecewise constant, with 
discontinuity points separated by at least e and jumps no larger than n'e (see Fig. [6] for an 
illustration). This implies that 


Vlv G SBV loc (R 2 -,R 2 ), 

consists of parallel straight infinite lines with normal separated at least by distance e, 
Dfi^ v — I C 2 + b ei/ 0 N v V} [j p , b ei , • = 0, b eiy G R 2 constant on each line, \b eu \ < n'e . 

(5.9) 


Let l o' — {X G uo : dist(X, duo) > L 2 e}. Then the plastic strain F^ is defined by 

FP U'= Dtp* L w /= I £ 2 U'+ J] b ei ® Nj H 1 • (5.10) 

j 

Here Jj are the finitely many segments which constitute the jump set of fi^ in uo, each of 
which has a constant jump b e j = [^e]. 

Note that we shall show below (Lemma |5.2| (ii)) that Dfi\ is uniquely defined in uo' and that 
this gives a unique definition of Fj? from fi e on O' \ ^(m+L 2 )e(X/c)^, where O' = {X G 

Li : dist(X,dft) > L 2 e}. 

(C3) For any dislocation X&, there is a map fi\, Lipschitz continuous and one-to-one, such that fi e 
can be expressed as 

<p e (X) = o <pf) (X) VX G B ce (X k ). (5.11) 

The function fi^ G S f HHi oc (R 2 ; R 2 ) is explicitly defined from the slip-line normals and the 
Burgers vectors. Precisely, in each core we are given N^k normals to the jump set Nj, 
c.f. (Cl), with slips b e j G R 2 , which obey |b e j| < n'e and b e j • Nj = 0 as in (C2). As a matter 
of convention, N 2 - will denote the counterclockwise 90° rotation of Nj, and the normals to 
the jump sets are chosen so that N j- points radially outward from the dislocation center, to 
the direction of the slip line, see Fig. |9| 

Given these parameters, we set 




Ndk n (SJ- Y A 

(X) = x - g Kj (X - x fc) ^ (x 2 ~ Xfc) , 


3 = 1 


where b*-(X) is the smoothed radial slip, 


b :, (x) = ^NVb e , min {iN, 1 } 


(5.12) 


(5.13) 
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Each function b*-(X) vanishes at the origin, is Lipschitz continuous and radial, in the sense 
that b*^(X) is parallel to X for all X. Further, on the part of the slip line outside B me it 
coincides with the “external” slip b e j, in the sense that 

= b € j for all t > me. (5.14) 

The angle Qj : M 2 —» [—27r, 2 tt] is a function which obeys 0j(te i) = 0 for all t > 0, 

X = |X|(cos0j(X), sin6f/(X)) for all XgK 2 , (5.15) 

and is continuous away from [0, oo)N^. Clearly 0j{tX) — 0j(X) for t > 0, and Qj is unique 
up to the value on the jump. 

We remark that cp^ is a SBV \ oc function, with jump set X& + Uj (0, oo)N^. Further, we recall 
that 

m = 2 N s ri, (5.16) 

c.f. (Cl). This conditions guarantees that cp^ is bijective. In particular, since all b*^ are radial 
we have <^(X& +Mv) C X& +Mv for all unit vectors v. An explicit computation shows that, 
for any unit vector v and t G M, 


vipZ-k + tv) = + V 


t — sgn(t) min < — 
1 me 


\ Ndk 

- 1 E N r 

} 3 =1 




6 jiy) 

2i r 


(5.17) 


This function is guaranteed to be bijective if (cp^ — X&) • v is a monotone function of £, or 
equivalently, d(cp^ — X&) • v/dt > 0. Since |b ej | < n'e, such condition will be satisfied by 
setting m as in ( |5.16 ). 

An explicit computation shows that 


|V^-I| <2. (5.18) 

Further, outside the inner core £> me (X&) the jump of (p^ has the same form as required in 
(C2), in the sense that 


Ndk 

D( Pel(R 2 \B m£ (X k )) = V^e'C 2 L(R 2 \B m£ (X fc ))+F be i 0 L(X fc +(me,oo)N^) • ( 5 - 19 ) 

3 =1 

(C4) The total deformation gradient and its elastic and plastic parts are respectively defined as 

F e = D(p e € M (ft; M 2x2 ) , (5.20) 

F® = Vip e E L 1 (ft; M 2x2 ) , (5.21) 

FP = I C 2 + b ei ® N i H 1 [ J: , nn , E M (ft; M 2x2 ) , (5.22) 

3 


where Fj? is set as I on and is defined as above. By (C2) and (C3), the jumpset J is 

a union of segments, i.e. J = U jjj, each with a constant vector b ej . The value of b e j inside 
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the core is taken for each j equal to that of the corresponding segment in the compatible 
region 7? ce (X&) \ Then Curl Fj? in fi / satisfies, c.f. Reina and Conti (2014:), 


Curl FUn'=E ±b ^°L(^ 


Jnfi'j 


(5.23) 


where the + sign corresponds to the endpoint of segment Jj pointed by N^, and the — sign 
to the opposite endpoint. The support of Curl Fj? |_q/ thus corresponds to the N e dislocation 
points. 

The fields F^ that result in a finite energy belong to the space BV (fl; R 2x2 ) and therefore the 
traces F^ + and F® - on each side of the jump set J are well-defined quantities by Theorem 
3.77 of Ambrosio et al. (2000). We then define F^ on the jump set J as 


F:(x) = F^Xl+irm vxej, 


(5.24) 


although the precise definition will not affect the limiting results. 


We first show that cpe is approximately Lipschitz. 

Lemma 5.1. Let pe : R 2 R 2 be a pure-slip deformation as defined in (C2), and set L — 
2 Ns (1 + n ') Ns . Then the following holds: 

(i) (fe is invertible, in the sense that there is p\~ l : R 2 R 2 such that p^~ l [p\ (X)) = 
<^(<^ _1 (X)) = X for all X G R 2 . 

(ii) For all X, Y G R 2 satisfying |X — Y| > e one has 

|^(X)-^(Y)|<L|X-Y| (5.25) 

and, similarly, for all X, Y € M 2 satisfying |X — Y| > e, 

\^ e ’-\X) - ^-\Y)\ < L\X-Y\. (5.26) 

(Hi) For any X G M 2 and any X € M 2 one has 

<^(£e(X)) C B Le (^ e (X)), (5.27) 

^~\B e (X)) C B Le (cpf~\X)). (5.28) 

(iv) Setting n = (1 + n') Ns one has, for all X in the jump set of p\, 

|[^]|(X)<ne. (5.29) 


Here and below we use the overhead tilde to denote quantities in the intermediate (compatible) 
configuration <^e(R 2 ). 
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Proof. Consider a single slip deformation <p^i, c.f. (C2). By construction, it can be written as 

^P ;1 (X) =X-^(X-N)N ± , (5.30) 

with <f e : R M piecewise constant, see Fig. [6| Further, the discontinuity points have a separation 
of at least e from each other, and the jumps of <f e are no larger than n'e. This implies | (j> e {s) — (j) e (t) \ < 
n'eff{J ( f )e n [s,£]) < n'(\t — s\ + e) and therefore 


|^ ?1 (X) - ^(Y)! + e < (1 + n')(|X - Y| + e) for all X, Y G R 2 . (5.31) 

The same holds for all Pe,w Taking the composition we obtain 

|(^P(X) - y>P(Y)| + e < (1 + ri ) Na (|X - Y| + e) for all X, Y G M 2 . (5.32) 


Setting L — 2 Ns (l + n') Ns concludes the proof of (5.25). Taking the limit Y —>► X + , X X - 
proves ( 5.29| ). 

To prove (i) we observe that (5.30) easily implies that <p v el is invertible, with inverse 

^;r 1 (X) = X + 0 e (X-N)N ± , (5.33) 


where <f> e is the same function as in (5.30). Then cp^ 1 is defined as the composition of the inverses 
of the single slips, and the estimate 


’ _1 (X) - <^p ,_ 1 (Y)| + e < (1 + ri ) Na (|X - Y| + e) for all X, Y G M 2 (5.34) 


follows. This concludes the proof of (i) and (ii) 
and (pT34h. 


Finally, assertion (iii) follows immediately from 

□ 


Now we prove that is well defined. 


Lemma 5.2. Let ip e , ip e e , cp 1 ^ and Fj? as in (C2-Cf). Then the following results hold: 

(i) The measure Fj? is uniquely defined from (p e . 

(ii) If Y is such that uo = B L 2 e (Y) C Ll \ U kB me (X.^), and cp e = ip\ o (p%, <p e = <p e e o (p f are two 
decompositions as in (C2) on the domain co, then on the smaller set B e ( Y) they coincide up to a 
translation, in the sense that there is c G M 2 such that <2?p(X) = <pP(X) — c for X G BJY) and 
^e(X) = *>®(X + c) for X G *£(B C (Y)). 


Proof. We first show that the second assertion implies the first. Recall that Fp = I outside 
O' — {X G O : dist (X,<90) > L 2 e}. Every Y G O' \ U^S( m+L 2 ) e (X^) has the property stated in 
(ii), therefore Fj? is unique on this set. Further, this shows that it is unique in the outer region 
£> ce (X&) \ B( jn j rL 2 ) e (Xfc) of each core, and therefore by (C3) also in the interior. This concludes the 
proof of (i). 

We now turn to the second assertion. Let Y be as stated. By the previous Lemma, p\ (B e (Y)) C 
B Le (y>P(Y)) and iff 1-1 ( B Le (^(Y))) C B L a e (Y), c.f. Fig. g We shall show that Dtp? = Dtp? on 
B e ( Y) by showing that J v and b e are uniquely defined from <p € on such domain. 

Uniqueness of By definition, cp e = <p\ o with ip\ G SBVio^M 2 ; M 2 ) and p\ Lipschitz 
continuous. Therefore, by the chain rule, c.f. Theorem 3.96 of Ambrosio et al. (2000), and since 
= I, 

D<f e = V<p e £ o C 2 + (p e £ (tp? + ) ~ &(<£-)) ® N H 1 L jr, , (5.35) 
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and J ( p e C p. For X G J7^p, <^e + (X) ^ ^ (X) and since p\ is one-to-one, then ^ ^<£e + (X)^ 7 ^ 

^e _ (X)^ . The function y> € is then not approximately continuous in X and Lt 1 p \ J = 0. 
Therefore, J7^p = and J7^p is uniquely defined from p e , up to null sets. 

Uniqueness of b e . We consider two decompositions p e = p\ o p 1 ^ = p e e o (pP in B L 2 € (Y). 
We write B L 2 e (Y) as the union of finitely many non overlapping connected open domains £>, 
i.e. (jj = U/3, such that B D J = 0. For each domain B , p^ and <pP can be written, by construction, 
as cpe (X) = X + cb and <£]?(X) = X + c#, with c# and c& constant vectors. Then, 


Ve(X) = (<pl o *£) (X) = {ip% O <p*) (X) = p\ (X + cb) = (fit (X + cb) ■ (5.36) 


Let d# = eg — eg. We proceed to show that all blocks B such that BDB e (Y) 0 have the same 
&&. For that purpose, consider two blocks B and B' such that they are in contact in the intermediate 


configuration defined by p^, i.e. (£> + eg) D (B f + C&) 7 ^ 0. Then, there exists X G £>, X' G S', 

such that X + c# = X' + c#/ and (X + eg) = ^(X' + c#/). Then, since <^ e = p\ o o 

^ (X + cb) = (p\ (X + cb) = ^ (x' + cbO = ^ (x' + C B ') • (5.37) 

Since p\ is injective, X + c# = X' + c#/. Subtracting this equation from X + c# = X' + c#/, 

one obtains that d# = d#/ or, equivalently, that c# — c#/ = eg — c#/ for every two pair of blocks 
in contact in Bl 6 (pP(Y)). Since its pullback to the reference configuration contains £> e (Y), one 
obtains b e = {(fej = l&i — b e on B e ( Y) and the result follows. □ 


The next Lemma shows that the decomposition of p% is unique, and will be used to prove 
det F p = 1. 


Lemma 5.3. 

L as in Lemma 


Let p e G X e as defined in (5.4). Consider Y such that dist (Y, dfi) > 2 Lre, with 
and B l 2 e (Y) satisfying B L 2 e (Y) n (U k B me (X k )) = 0. If <p\ and dp, 


5.1 


P are two 


deformations on B L 2 e {Y), satisfying Fj? = Dp\ — Dp\, each with its decomposition of the form 


( P€=V%,N. 0 - 0 V > e,v 0 - 0< f%,l’ 

$% = <f%,N s 0 - 0 ° - ° $*,1’ 

then, there are constant vectors c v , v — 0, ..,N S such that 

p>l u (X) = <piAX- c„-i) + c„, 

D Vlu (<Pe, v -l ° ••• ° <l) = D( P% ($£„-1 ° - ° vil) 


(5.38) 


(5.39) 


in B e ( Y). 

Proof. Consider two possible plastic deformation mappings p^, p^ on B L 2 € (Y) and their decompo¬ 
sition as in (5.38). By Lemma 5.2, Fj? = Dp\ — Dp\ on B L 2 e (Y), and by integration one obtains 
Pe = pP + cw s , where cn s corresponds to a rigid translation. Define p v eN ,(X) = £V( X ) + c Ns . 
Then, 

vIn s °vIn s -i° ■■■°vIi = vIn s °V% s -i° on B L 2 e (Y). (5.40) 

The last deformation applied in the sequence corresponds to slip along jump sets with normal 
N n s . Since the normals N^, v — 1,...,7V S are all distinct, the jump set with normal N n s on 
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Figure 7: Choice of domains to prove the uniqueness of in B e ( Y), for B L 2 e (Y) C D\U &£> me (X&). 


configuration j o ... o ( B L 2 e (Y )) will necessarily traverse regions where F p = I. The 

segments in those regions will suffer a translation upon pullback to the reference configuration, 
and will thus have that same normal Ny, each with a given Burgers vector. The ensemble of 
segments with normal Nw s will, by construction, form a set of straight lines in the deformed 
configuration with a constant Burgers vector per line, potentially different between lines. Thus, 
D <Ns is uniquely characterized from F p on the set < (B L 2 e (Y)), which implies 1 = d vIn] 

and 


<W!(X) = &£(*) + C N S - 1 = <%£(X - C Ns ) + Cjy g —i 
^e,JV s ( X ) = ^e,JV s ( X ~C Ns -l) = <pf,N s ( X ~ c N s -l)+ C Ns 


on D Ns , 
on ^n!( d n s ), 


(5.41) 

(5.42) 


where Dn s = F>Le(^e(Y)) C cfe ( B L 2 e (Y )) and cw s -i and cn s are constant vectors. We define 
E Ns = (fe’- 1 ( D Ns ) so that 


^e,AT s -l ° - ° ‘fill = £1^-1 ° - ° £e,l + c iV s -l on ^JV. • 


(5.43) 


Defining at s -:l(X) = ^ iVs _ 1 (X) + c/v s -i, it is readily obtained 

• o^e,i o n£, iV s - (5-44) 

The above procedure may be iterated for v — N s — 1 successively till v — 1, obtaining for all z/ 


<JV S -1 ° <JV S - 2 ° 


° <1 = ° ^ 6 ,iV s -2 ° 


^(X) = ^(X) + C„ 

<,(X) = <^(X - c„_i) = £P„(X - cw) + c v on ^(A,), 
H^(X) = D<pl„(X - c,,-!) = Dft tV (X - c v _i) on ^(A,), 

<*,-1 ° ‘fills—2 ° ••• ° <1 = £e,„-l ° ‘fills -2 ° ••• ° £e,l> ° n 


where £)„ = B 


Le/( 2 + 2 n') N s~ 


<fie 


-°<i)(Y) ),E„ = (<p P 






(5.45) 

(5.46) 

(5.47) 

(5.48) 

and c„ are con¬ 


stant vectors. By (5.31), D u C ^^(-A+i)- Note that for v — 2, the last expression gives 
VjP-^X) = y^yx), and therefore Co = 0. 

The composition of two consecutive deformation mappings then gives 


‘Plis(‘filis-i{ x )) = ‘filu (<„-i(X) - c„_i) +c u = Cp\ v ($£„_! (X - C^_ 2 )) + C, 


(5.49) 
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and the composition of the first v mappings yields 


=¥ , £,, 0 ^-i 0 -°¥’e,i( X_c o) + c ‘' °n Ex, (5.50) 


where cq = 0. Thus, from (5.47), 


D( Pt,v (v>e,u-l ° - ° ^6,1) = D( P% (<^-l ° - ° <1 - C ^-l) = D( P% (fie,v-l ° - ° &e,i 


We note that B e (Y) C Fi, and E v C E^+i, and thus Eqs. (5.39) hold for all v. 
Lemma 5.4. Let J be the jump set of cp e E X e . Then, the following results hold: 

(i) There exists A > 0 such that for any ball £> r (X) C M 2 of radius r > e 

Ar 2 

\jnB r (x)\ < —. 

(ii) There exists C* > 0 depending only on such that 


on E\. 
(5.51) 

□ 


(5.52) 


(5.53) 


We recall that by definition J C(l. 


Proof Consider first a ball B e (Y) such that £>2Le (Y) H (U&i? me (Xfc)) — 0. In this ball, the 
representation of from (C2) holds. Since the are piecewise translations and the normals 
N v are all different, the jump set of is the union of the counterimages of the jump sets of the 
individual p 




u^r 1 ^ 


<2 


u 


U(w 


el _1 



For every v, let Y v = (p^v-i o • • • o y>^ 1 )(Y). Since the jump set of pe,v consists of parallel lines 
with distance at least e, we have 


PB Le (Y u )\ < (2L + l)2Le 

and, recalling that {ip\ v _ x o • • • o ^p v e l ){B e ( Y)) C B Le ( Y„), 

l(^r 1 o---o^-_ 1 i)(^pJnS e (Y)| < (2L + l)2Le. 

Therefore |J^p n S C (Y)| < N S {2L + l)2Le. 

Consider now a generic ball £> r (X). We cover B r (X) \ (U&£> ce (X&)) with balls uoi = B e (Yi) with 
finite overlap such that £> 2 Le(Y/) H (U kB me (X.A) = 0. The number of those sets ui is bounded by 
Air 2 /e 2 for some constant A\. Similarly, by (5.6) the number of cores intersecting B r is bounded 
by A 2 r 2 /e 2 for some constant A 2 , and the total length of the jump set in a core is not larger than 
ceN s . Then, 


\J n Br\ < ^2 \x n coil + J 2 \x n (B r n B ce (x k )) | 

l k 

2 2 2 
< A\— y (2 L + 1) 2 Le + A 2 ^N s ce < A —, 


(5.54) 
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which proves (i) 

For assertion (ii), since ft is a bounded domain, it suffices to cover it with balls of the two types 
used in (i), or to choose in (i) a ball which covers ft. □ 


The next Lemma will be used to show that the elastic energy associated to cores vanishes in 
the continuum limit e 0 . 


Lemma 5.5. Let cpf be as defined in (5.12), and N e and S ce (X^) as defined in (Cl) and (C3). 
Then there exists C c > 0 such that 


N e r 

Y / IV^I 2 dX < C c e. 

k jB ^ k ) 


(5.55) 


Proof. From (5.18) and scaling (5.5), the result is readily obtained 

N e „ N e 


Y t \V<Pe\ 2 dX < 16Y\B ce (X k )\ < C c e 
k 7 - 


(5.56) 

□ 


In closing this Section we show that the decomposition we assumed for the plastic strain is 
generically admissible. This is not used in the current argument, but illustrates the generality of 
our assumptions. 

Lemma 5.6. Consider an arbitrary F p G C foo (fl;R 2x2 ) such that detF p = 1, Ff x 0 and 
Ff 2 7 ^ Ff ± . Then F p admits a unique description as the composition of the following three sequential 
simple shear deformations 


F p = (I - 7 e 2 <g> ei) (I + p ei 0 e 2 ) (l + /a e : 
with e 3 = y/2/2 (ei + e 2 ). In particular, 

F* - 1 


fi — 2 


11 

FI2- 




r 1 = l + F? 2 -F» 1 

1 + F 2 p i-F 2 p 2 


7 = 


F?2 ~ 


where Ff ± , Ff 2 , F 21 and F 22 are the components of F p . 


Proof. Expanding the product in (|5.57|) one obtains 
F p = 




F p 

r 21 


F?2 


F P 

r 22 


+ 7f 


1 e 3 


-f + 7 + 7§ 


-7 + 7f “ 2 77 +2 7 2 “ 77 - 77 2 + 1 + 2 


(5.57) 


(5.58) 

(5.59) 

(5.60) 


From the expressions of F)) = 1 + ()(r/ — 1 ) and Ff 2 = 7 + ^(7 — 1), one obtains 77 = 1 + Ff 2 — F[\■ 
If Ff.' 2 / Fj :> |, then, plugging the result into the equation for Fj) we get Ff) = 1 + ^ (Ff 2 
and so (5.58). Finally, the expression for F 2 \ simplihes to 


Ffv 


7, = f +7[-f(»-l)-l]=f-TFT,. 


(5.61) 
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Using the previous result for /i, and the fact that detF p = 1, one obtains 


iK = 


7 ?P _ 7 ?P 7 ?P 7 ?P 

^11 ^ 11^99 + ^ 11 ^ 9-1 


lP 22 


IP 21 


f & 


which simplifies to (5.60) if Ff x ^ 0 


(5.62) 


□ 


We note that the choice of the orientation of the three slip systems was arbitrary and different 
orientations may be chosen. The conditions for obtaining a unique representation under such choice 
will then differ from that of this lemma. 

6 Main results. 


Proposition 6.1. Let cp e G X e 


as defined in (5.4), F\ 


and Fg as defined in (5.21) and (5.22) 
respectively and L as defined in Lemma 5.1. Assume sup e E e (cp e ) < oo. Then, there is C\ depending 
on sup e E e ((p e ), such that for any set K compactly contained in and for any e G (0,1) satisfying 
e < dist(dfl,AT)/( 2 L 2 ) 

\D<Pe ~ F®FP|(if) < Cie 1 / 2 . (6.1) 

Proof. We decompose the proof into three steps that concern, respectively, the dislocation cores, 
the compatible subdomains away from the dislocations and, finally, the total domain. These in¬ 
termediate results will often use a combination of the Boundary Trace theorem (Theorem 3.87 
of Ambrosio et al. (2000)) and the Poincare inequality (Theorem 3.44 of Ambrosio et al. ( | 2000 | )) 
for BV functions, that we make precise here. Note that the constant in the Boundary Trace 
theorem depends on the domain. If the domain is a ball B\ of radius 1, the theorem asserts 
SdB! I f\dV} < Cn||/||.BV(.Bi) = Cti f Bi \f\dX + C t i\Df\(Bi) for some constant C t \ > 0. If B x is 
replaced by £> e , where B e represents a ball of radius e, a scaling argument yields the inequality 
IdB \f\ dnl — Ib \f\ dX + Cti\Df\(B e ) for the same constant Cti as before. In particular, if 
/ GW(B e ), then 


L 


c tl 


\f-f\dH 1 <C tl \Df\(B e ) + ^ f |/ 

dBe e JB e 


f\dX < C tp \Df\(B e 


(6.2) 


where / = \Bf\~ 1 f B fdX and the second inequality follows from Poincare’s inequality, see also 
Ambrosio et al. ( 2000 ) [Remark 3.45] for a corresponding scaling argument. The constants Cti and 
Ct p are independent of e and of course, independent of /. The above result is also true for J D B e , 


/ 

JjnB e 


1/ - JldU 1 <C tp ,\Df\(B e 


(6.3) 


as long as J PB e is the union of at most M segments, where the constant C tp ' depends only on M. 
This estimate holds for / on the jump set being defined as either one of the two one-sided traces / + 
and /“, or as the average (/ + + /“)/ 2 . To see this, let C p \ be the constant entering the estimate 


corresponding to ( 6 . 2 ) for the half ball, which does not depend on the orientation of the half ball. 


Given a segment 7 C £> e , consider the diameter 71 parallel to 7, and let 7' be the projection of 7 
onto 71 . By the reverse triangle inequality and the fundamental theorem of calculus, working for 
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Figure 8 : Sketch of the strategy to bound the absolute value of a function over a segment in a ball, 
with the variation of the function over the ball. 


notational simplicity in coordinates such that 7 is parallel to ei, 

[\ / - f\dn l - [ 1/ - ,/W = [ (1/ - f\(X,Y 2 ) - 1/ - /|(x Yi)) dft 1 

A Jy J\x i,Xoi 


< f (|/-/|(X,Y 2 )-|/-/|(X,Y 1 )) 
•/[*!,* 2 ] 


dft 1 < 


7 1 / 

•/[X 1( X 2 ] d[y!,r 2 ] 


D 2 (f - J)dH ] 


dn 1 < [ D 2 (f - f) 

J Dry 


[Xi,x 2 ] 
dX 


((/ - /)(X,Y 2 ) - (/ - /)(V U)) 


(6.4) 


dl-L 1 


where D 2 denotes the distributional derivative in direction orthogonal to segment 7 , and (X\, Y 2 ), 
(X 2 , Y 2 ) are the endpoints of 7 , (Xi, Yi), (X 2 , Yi) the endpoints of 7 ', see Fig. |8j Therefore, 

[ \f-]\dU l < [ |/-/|d'W 1 + |A/|(Z) 7 ) < [ \f-j\dn l + \Df\{B e )<{C pl + l)\Df\{B e ), (6.5) 

so that for M segments, (6.3) holds with C tp > = M{C p i + 1 ). 

(SI) In this first step we show 

\D<p e | ^JfUXfc)) <077 ( 6 . 6 ) 

|FfFP| (\j£ C£ (X fc )j < C 3 e 1/2 . (6.7) 


We start with ( 6 . 6 ). The displacement jump at X G Jj H B ce (X.jr) is 

b £ (X)] = y>+(X) - y>r(X) = 7+(X) - ¥>+ (Xfc)) - 77(X) - <7 (Xfc)) ( 6 . 8 ) 


where we have used <pf (X&) = <£>“ (X&), c.f. (C3) in ( |5.4| ). We further note that between X& 
and X either above or below Jj we do not cross any other jump set, and therefore Dip e = F 
on both sides. Since has a trace on each side of the jump set, we obtain 


¥>t(X)-¥>+(X fc ) = /^Ff+dX', 
dx fc 


(6.9) 
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B ce (X k ) 



Figure 9: Example of a dislocation core £> ce (X&). The dislocation is centered at X& and is the 
boundary of a finite number of slip lines Jj, 1 < j < N^k < N s . 


and the same on the other side. Inserting in (6.8), 


vt (X) - & (X) = /°W - FD dX !, 

■/Xfc 


( 6 . 10 ) 


which implies |[<^ € (X)]| < \DF^\(B ce (X.k )). Integrating over one of the slip lines inside the 
core, which have length ce, and then summing over all of the slip lines, we obtain 


Ndk n 

X / M dH 1 < N s ce\DFt\ (B ce (X k )), 

-• JJjnB C€ (x k ) 


( 6 . 11 ) 


where we have used N^k < N s , c.f. (Cl) in (5.4). Summing now over all dislocation cores, we 
obtain 

N f N dh „ / N e \ 

I foil dH 1 < N s ce\DF e e \ \jB ce (X k ) < C±eE e . (6.12) 


N dk 

EE 


k j 


ST j nB ce (Xfc) 


We turn to the absolutely continuous part of the gradient of cp e , which by Holder’s inequality 
can be estimated as 


Ne N e / N e r 

X / l V ^l = X / l F eI < V^ceivy 2 X / 

, 4 B ce (Xfc) , is ce (X fc ) \ k Je 


k ” B ce (Xfc) 


|F®| 2 dxj 


1/2 


(6.13) 


Recalling the quadratic growth of the energy, c.f. ( |5.3|), and the definition (5.2) we see that 
the parenthesis is bounded. Since < C/e (recall ( pn5| )), this proves ( |6.6[ ). 

We now turn to (6.7). By (5.24) and the bound \b e j\ < n'e < ne, c.f. (C2) and (C3) in (5.4), 

r Ndk r 

\F e e F p e \(B ce (X k ))< \F e £ \dX + neS2 |F e e \dH l . (6.14) 

JB ce (X k ) - 4j;-nB ce (Xfc) 
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In the second term we use (6.3), so that it becomes 


N dk r 


|F®| dU 


N dk r 

L 


j J J^n-B ce (X/ e ) j j nB ce (Xfc) 

where F e denotes the average of F® over £> ce (X&), which obeys 

1 


|F*| dU 1 + C 5 ne\DF e £ \(B ce (X k )) , (6.15) 


|Ff(X*)| = 


7rc 2 e 2 


[ F® dX < f |F: 

JB ce (x fc ) ttc e J Bce (x k ) 


®| dX. 


Therefore 


77 TV f C 

|F*FP| (B ce (X k )) < (1 +-f) / |F®| 7X + C 5 ne|HF®|(S ce (X fc )). 

WC JB ce (X k ) 


(6.16) 


(6.17) 


After summing over k, the first term is bounded as in (6.13), the second as in ( 6 . 12 ). This 
concludes the proof of (6.7). 

N e 

(S2) Next, we cover the domain K \ (J £> ce (X&) with finitely many balls £> e (Y/), with Y i E K, 

k 

such that the larger balls B L 2 e (Yi) do not intersect the dislocations cores B me (X&) and have 
finite overlap. Such a covering is possible in view of Eq. (5.6) and, by the definition of K, 
B l 2 e (Y / )cO , ca 

We shall now show that for any l one has 


I Dtp, - FfF*| (B e ( Y,)) < C,e\D¥l\(B he (Yi )), 


(6.18) 


N e 


with h— 1 + Ln < L 2 . Summing over all the balls will give the estimate on K \ (J £> ce (X&). 

k 

Fix l and let and (fe be as in (C2), applied to the set £> 2 L 2 e (Y/). Then, p e = o 
Taking the gradient we obtain Y e e — \7p e e op^ and 


D<p e = F \C? + YA<P J » N ? U 1 [j. in B ke ( Y; 


(6.19) 


Therefore, recalling (5.10), 


l^ e -F«FP|(B e (Y / ))= [ 

J E 


BeiY^nj 


\l<p e }-F e e b e \dH l , 


( 6 . 20 ) 


where b, takes the constant value b CJ on each segment J :j C B e (Yi) fl J . For X 6 B c {Yi)C\Jj 
we write X+ = <^? + (X), X = ^“(X) and X+ - X“ = b ej . Then, 


bJ(X) = ^(X+)-^(X-)= B V<p e e dX= f 1 F e £ (^- 1 (X~+tb ej ))b ej dt . 

JX- JO 


( 6 . 21 ) 


The value of Ff on the jump set was d efine d in ( 5.24 ). Since S — [X + ,X ] is a segment 
of length no larger than ne, by Lemma 5.1, pe’ L (S) C Bh e (Yi). Furthermore, 
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consists of the union of at most M = (1 + n) Ns segments. Indeed, each ip 1^7 1 is piecewise a 
translation, and a segment of length no larger than ne crosses at most |_^J + 1 discontinuity 
points, hence gets decomposed into at most |_^J + 1 segments, each of length no larger than 
ne (plus possibly a finite number of points, which may be ignored). Iterating at most N s 
times gives the result. Then, by (6.3) one obtains 




'vr \s) 

where F e is the average of F® over S^ e (Y/), 

Fl-^- 


|F® - F*| dU l < C 7 \DF e e \(BteiYt)), 


\Bhe(yi)\ J Bhe(Yl ) 


F tdx. 


Therefore (6.21) and (6.22) yield 


bJ(X) - FXI < [ |F« - f:| dU 1 < C 7 \DF e e \(B he ( Y,)) 


( 6 . 22 ) 


(6.23) 


(6.24) 


We write (6.20) as 


p¥» e -®1FP|(B e (Y,))< [ 

Je 


Be(Y t )nj 


I^J-F^bel^ + lbel [ \Ft-Fl\dU 1 . (6.25) 

JB e (Y,)nj 


Since, by (6.3), 


lB e (Y t )nj 


iFl-Flldn^CslDFtKBUYt)), 


(6.26) 


this concludes the proof of (6.18). We remark that by (C2) and the proof of Lemma 5.4 
\B e (Yk) fl J\ consists of a finite union of segments bounded independently of e. 

(S3) We combine the results for the dislocation cores in (SI) and each B e {Y{) in (S 2 ), which, by 
construction, have a finite overlap. Then, 


\DtPe ~ FgFP| (K) < Oe 1 / 2 . 


(6.27) 


□ 

Theorem 6.2. Let p e G X e , F e ,F® and F^ as defined in Section^ and sup e E e (p e ) < oo. Then 
there exist F p , F e , ip and a subsequence (not relabeled) such that 


j?p A pp 

in M (fl;M 2x2 ) , 

(6.28) 

pe pe 

in L 2 (Q-R 2x2 ), 

(6.29) 

pe pe 

in BV (fl; R 2x2 ) , 

(6.30) 

pe pe 

in L 1 (ft;M 2x2 ), 

(6.31) 

p e -v e ^p 

in L\ oc (f2;M 2 ) , 

(6.32) 

Be = Dp e A 

Dp = F in Mioc (fi; M 2x2 ) , 

(6.33) 

Curl F p Ln'- 

- Curl F p in M (fi; M 2 ) 

(6.34) 
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as e ^ 0 , where v e is a sequence of rigid translations and ft' = {XE(1: dist (X, an) > L 2 e}. The 
support of the measure CurlFg |_n' thus exclusively consists of dislocation points. Furthermore cp is 
approximately continuous TL 1 -almost everywhere. 


We recall that fi e is said to weakly* converge to fi in A / 1 i oc (^) 5 if 

lim / fdfi e = / fd/j, 


6 —t*0 


(6.35) 


for all / G C'c(fi), c.f. Definition 1.58 of Ambrosio et al. (2000). 


Proof. Convergence of Fj?. 

By the definitions of F^ and J , and the bound on b e j ( 5.29| ), |F|? |(D) < \/2|^| + |b e j|| < 
\/2|D| + ne\J\. Therefore sup e |Fe |(D) < oo and (6.28) follows by weak* compactness (Theorem 
1.59 of Ambrosio et al. (2000)). 

Convergence of F|. 

By the growth condition of the elastic energy (5.3), sup e ||F||| L 2 (^) < oo. Then, (6.29) follows 
by weak compactness (Theorem 1.36 of Ambrosio et al. (2000)). Furthermore, since D is bounded 
and the energy imposes sup e |DF^|(D) < oo, we obtain sup e ||Fg Ws vin) < oo. Hence (6.30) and 
(6.31) follow by BV compactness (Theorem 3.23 of Ambrosio et al. (2000)). 

Convergence of cp e and Dcp e . 

Consider a fixed closed ball B compactly contained in D, and define v e = j^j f B F> e dX. Next, 
choose K' CC D and define K = K'UB. Then for e < dist (AT, dQ) / (2 L 2 ), we show that \D<p e \ (K) < 
oo by decomposing the domain into the dislocation cores and the remainder of the domain, c.f. proof 
of Proposition 6.1 By (6.6), 


'N e 


I^V.I lK(x») < 


(6.36) 


Next, we cover K \ (J^ B cc (X*,) with balls of radius e, B e (Yi), l € N with finite overlap, as in 
Proposition 6.1 (S3) and we further define 


nm = -5 


F; dX. 


7re^ 


lB e ( Y,) 


(6.37) 


By Proposition 6.1, \Dip e — F®Fe | (K) < C\e l / 2 . The bound on |FfFe|(/T) is obtained with a 
similar strategy as in Proposition 6.1 For each ball, since h C] \ < ne we have 


|F®FP| (B e (Y t )) < [ |F°| dX + nef |F HdH 1 , 


(6.38) 


average via (6.16) as in (6.17) gives 


and J fl B e (Yi) consists of finitely many segments. Using (6.3) as in (6.26), and estimating the 

(6.39) 


|FfFP| (B e (Yi)) <C 9 [ |F®| dX + C 10 ne|Z)F®|(B e (Y/)). 

JB t { Y,) 

Summing over all balls, 

|F«FP| (K \ (U k B ce (X k ))) < C u f |F«| dX + C 12 e\DF e e \(n) 

J Q 


(6.40) 
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and therefore, recalling ( |6.36| ), 

\D<Pe\ (K) < \D<p e \ (u fc B ce(Xfe) ) + | Dip, - F°F£| (K) + |F«F£| (. K \ U k B ce ^ k) ) 

< C 2 e 1/2 + C x e l/2 + C n f |F®| dX + eC 12 \DF e e \ (fi). 

J Q 

Therefore sup e \Dcp e \ (K) < oo. Then, by the Poincare theorem for BV functions, c.f. Theorem 


(6.41) 


Ambrosio et al. 


( 2000 ), sup e \\ip e — V^elUhiT) < 00 5 where cp e is the average of p> e 


over 


3.44 of 

K. Since B C K it follows that |v € — <p e \ is bounded as e 0. Then, as a result, by BV - 
compactness, c.f. Theorem 3.23 of Ambrosio et al. (2000)) and the definition of weak* convergence 
in BV (Definition 3.11 of Ambrosio et al. ( 2000| )), ( |6.32| ) and (6.33) follow. 

Continuity of the limiting solution 

From the previous results, cp e — v e —» cp in L 1 1 oc (D). Then, by lower semicontinuity of the 
variation measure, c.f. page 172 of Evans and Gariepy| (|l991), 


\D(p\(B r ) < liminf \D(p e \(B r ) \/B r CCfi, 

e —^0 


(6.42) 


where B r denotes a ball of radius r. We write £>2r for the ball with the same center and radius 2r, 
and assume £> 2 r CCfi. 

Let us consider a set K CC D, such that £>2r CC iF, and take e < dist(<90, K)/(2L 2 ). Then, 
defining K e — K\ (l4£> ce (X&)), we obtain by Proposition 6.1, (6.36) and the definition of F? 

/ N e \ 


\Dip e \(B r ) < \Dcp e \(B r f) K e ) + \Dcp e \ B r Pt[jB ce (X k 


N e 


(6.43) 


< I D<Pe ~ F e eF p e\(Br n IQ) + |FfFP|(B r n IQ) + \D<p e \ (£ r n |J J B ce (X fc )J 

< C^e 1 / 2 + \F e e FP\(B r n K e ). 

To estimate the last term we cover B r nK e with balls of radius e which do not intersect the cores, as 
done in (6.40). The balls have finite overlap and are all contained in the larger ball £>2r• Therefore 

\D<p e \(B r ) < QQ' 2 + C 14 ||F^|| il(B2rnXe) + C lb e\DF e e \(BQ . (6.44) 

a 

(6.45) 

(6.46) 

(6.47) 


Passing to the limit, since (6.31) implies ||F||| L i( B2r ) —* ||F e || i i( B2r ), we obtain 


\DQ(B r ) < lirnrnf \D<p € \{B r ) < C 14 ||F e || Ll(B2r) < C 14 20Fr||F e || L2(S2r) 


Therefore 


Since, by (6.29), F e € L 2 (ft), 


\DQ(B r ) < C 16 r||F e || L2(S2j . ) . 
\DQ(B r 


lim inf 

r—>• o 


= 0 . 


As a result, c.f. Proposition 3.92 of Ambrosio et al. (2000), the jump set of Dp> has T~L 1 measure 
zero, and (p is approximately continuous / H 1 -almost everywhere. We note that it may contain a 
Cantor part. 
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Convergence of Curl F^ 


By the scaling (5.5) and the bound on the Burgers vector (5.29), Curl F|? |_o' is equibounded. 
Therefore, by weak* compactness (Theorem 1.59 of Ambrosio et al. (2000)), 


CurlFP[o' A G inM(ft;R 2 ). 


(6.48) 


Since Fj? also converges in Al, we have G = Curl F p . Indeed, for every G C^ft), one can 
choose e such that supp (j) CC ft'. Then, c.f Appendix A, 

lim [ fid (Curl F p - Curl F p ) = lim [ fi d (Curl F p - Curl F p ) = lim - [ Dfixd( F p - F p ) = 0, 
e ^° Jfl JQ' JQ' 

(6.49) 

where the cross product is considered between Dfi and each column of (F|? — F p ). 

□ 


as m 


(6.50) 


Theorem 6.3. Let (p e , F e = Dcp e , F| ; F£ as defined in Sectionand F = Dip, F e and F p 
Theorem 6.2. Further, sup e E e (p e ) < oo. Then, 

F F e F p . 

Additionally, F p G L°° (ft; R 2x2 ), F G L 2 (ft;R 2x2 ) and p G IT 1 ’ 2 (ft; 

Proof. In the following, we show Dp e F e F p in the sense of distributions. By Theorem 
uniqueness of the limit, this implies (6.50). 

We begin by defining fi e as the mollification of F^ with a kernel rj e G C^°(B e ) satisfying 
f B dedX = 1 and rj e (X) = rj e (-X), i.e. 


6.2 


and 


^e(X) = (F^ * r ?£ )(X) = [ rj e (X - Y) dF p (Y). 

J B e 


(6.51) 


Then, for any K CC ft, and e < dist(dft, K)/(2L 2 ), one has, by r] e < Cnj^XBe, (5.52) and 
( 5.29| ), that || / jL € II l°° (K) is equibounded, the bound does not depend on K. Therefore there exist 
G L°°(ft; R 2x2 ) and a subsequence satisfying \i e \i in L°°(K; R 2x2 ) for all K CC ft. Next, we 
show that fi — F p . Let fi G C£°(ft). By choosing K such that supp(</>) CC K CC ft, and choosing 
e < dist(dft,AT)/(2L 2 ), 


[ <MdF p -d/x e )= [ 0(dF p -d/x e ) = [ fi(dF^ e -d(F^rj e )) = [ (fi-fi*rje) 
Jn JK JK JK 


dF p 


< sup | (f) — (f) * 77 e ||F p |(ft) —y 0, as e —>► 0, 


(6.52) 


since |F^| (ft) equibounded, c.f. proof of Theorem 6.2, and fi*r) e fi uniformly as e 0. Therefore 
!i = F p and F p € L°° (Q: M 2x2 ). Furthermore, writing for simplicity, with abuse of notation, 
Dip e dX and Fl’dX for dD<p e and dF p , 


[ <f> (D(p e - F e F p ) dX = [ <f> (D(p e - F e F p ) dX 

in Jk 

<[ \4(Dcp e -F e e FV)\dX+ [ 0(FfF p -F ^ e ) dX 
Jk Jk 

+ [ \ct>{Fhx e -F e ix e )\dX+ f (j) (F e /u, e — F e F p ) dX 
Jk Jk 


(6.53) 
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The first term, f K \<f> (Dcp e — F|Fe) | dX, tends to zero by Proposition 6.1 The third term, 

Ik |^(F°/* e -F e /* € ) \dX, tends to zero by the strong L 1 (K ) convergence of F®, c.f. ( |6.3lD , and 

fi e G L°°{K). The last term, f K (j) (F e /x e — F e F p ) dX , tends to zero since </>F e G L l (K) and fi e 


weakly* converges to F p in Finally, in the second term, f K <f> (F|Fe — F|/x € ) dX 


, we add 


and subtract the identity, use I * r\ e = I and a change of variables in the second term, to get 


[ 0(F 

JK 


F p -F MdX = 


L 


Ik 


<t>K (FP-I)-0F«((FP-I)*7, e ) 
^Ff(FP-I)-((^)*,7 e )(FP-r 


dX 

dX 


(6.54) 


0F® — (</>F®) * rj e 


'N dK 


JKnj L 

We cover the jump set J with balls B e (Yi) such that the double balls I? 2 <:(V) are contained in Q 
and have finite overlap. We use (6.3) in each ball, both for (j) F® and for (</>F®) * rj e , the chain rule, 
the quadratic growth of the elastic energy and the energy bound, 


f ^(FfFP-F l» e )dX <J>e f 

JK i JB 




B e (Y,)nj 

<Y, neC is\ Z?(«|(S 2e (Y { )) 


dn 1 


(6.55) 


< C 19 ne ^|| L oo (n) |£>F||(fi) + \\Dct> || L =o (n) j 


F?| dX < C 20 e. 


Therefore Dp e converges distributionally to F e F p . By uniqueness of the limit, then F = F e F p . 
Further, since F e G L 2 (fl; R 2x2 ), and F p G L°° (Vt\ R 2x2 ), it immediately follows that F G 
L 2 (fl;R 2x2 ). By the Poincare inequality, (p G L 2 (fl;R 2 ). 

□ 


Theorem 6.4. Let F p be as in Theorem \6.ty sup e E e (<p e ) < oo. Then detF p = 1. 


Next, we cover fl 2 \ U&£> ce (Xfc) with balls B e (Yi ), such that Y ] G fl 2 an d |Y/ — X&| > ce for all Z, k. 
Note that this cover is possible by the constraint on the minimum distance between dislocations, 
c.f. Eq. ( |5.6| ), and our choice of fl 2 . We fix a mollification kernel rj a G where = e/10. In 

£> 2L 2 e (Y/) we choose a decomposition as in (C2), ° ••• ° Vci- We define, in the smaller 

ball S e (Y/), the regularized plastic deformation by 


Proof. Consider the set 0 2 = {X G : dist(X, dfl) > 2L 2 e}, with L as defined in Lemma 


5.1 


<p p = 
t ea 


(^Ins * w) ° ••• ° * Va) o ... O ^Pj * Tfc) 


(6.56) 


We observe that is smooth, and det I? (<£>?,i/*T7 a ) = 1 implies det = 1 everywhere. Further, 

since the jumps of each are bounded by n'e and separated by at least e, c.f. (C2), we obtain 
\<Pe,v * t\ol ~ ( Pe,i' | < n ' e and \D (<Pe,v * Va) | < cn'e/a. Iterating we deduce that \(fiea ~ Pel < C2ie 
and \D(p% a \ < C 22 . 

We now show that Dp v eOL is uniquely defined, independently on the choice of l and the decomposi¬ 


tion of cpe . Let <£ p = V^o-o^ibea different decomposition. By Lemma 5.3, there are constants 
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Figure 10: Cover of Ft" with balls B eP centered at the dislocation points X& (grey balls), and balls 
B eP /2 with finite overlap (blue balls) that do not overlap with the dislocation cores B me (Xk). 


c v such that yjj ^X) = (^ (X-c^-i)-^, which implies ((p%*r) a )(X) = (^ l/ *7/ a )(X-c i/ _i) + c l/ . 
Iterating as in ( 5.49||5.50 ) gives Dcp p a — D<p p a . We define F^ a as Dcp^a on UiB e (Yi) and zero else¬ 
where. 

Next, we choose a constant 1/2 < p < 1, and define FI" = {X E Ft : dist(X,<90) > e p }. We 
now construct a cover of Ft" with balls of radius of order e p . Precisely, we first select balls B eP {X.^) 
centered at the dislocation points, and then finitely many Y/ E H = {X E Ft : dist(X, dFt) > 
e p /2} \ UkB eP / 2 0^k) suc h that H C U/£> e p/ 4 (Y/) (this is possible since H is compact). Let (j) f k E 
C^°(S e p/ 2 (Xfe)), (j)i E C£°(i? € p/ 4 (Y/)) be a partition of unity on {X E Ft : dist(X, dFt) > e p /2}. We 
define Q' k = (j) r k * rj eP / 4 and 9 l = fa * r] eP / A . Then J2k^k + ft = 1 on ^ G C^°(S e p(X fc )), 
ft G C™(B eP/2 ( Y Z )), and |£>^|, |£>ft| < 

For (j) E C£°(fi), consider iF such that supp(</>) CC iF CC Ft and take e < (dist(iF, dFl))^^ p \ so 
that K C Ft" 

as above F^dX for dFj?, and uji — £> e p/ 2 (Y/ 


Then, using the previously obtained relation ft + ^ = 1 one obtains, writing 


lim 

6—^0 


[ (F» a -¥*)</> dX — lim [ (FP q -FP )<j>dX 
Jn e ^° Jk 


< I™ X I [ 0i(F?* - F^dX + lim V I [ 0' k (F% 

H0 , 


FP)0dX 
(6.57) 


In the first term we integrate by parts and use the previously obtained relations |<pP a — ipf | < C 2 1 e 
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Figure 11: The figure shows the support of F 9 in red (left image) and F 9 a (right image) near a 
given dislocation core £> ce (X&). The support of CurlF^ is X&, whereas the support of CurlF^ will 
be Ndk < N s diffuse dislocations located at the boundary of the domains B e in which a global 
definition of F 9 a was obtained. 


and \D0i\ < 9 , 

‘ToEl jpds. 

i Jaj ‘ 


FP )(j>dX 


= Jim^ / d^D^-Dcp^dX 

£ i Ju i 

= XI I [ - <Pe) D (0l<i>) dX 

< lim C*24£ (e + H-D^Hl 00 ) = 0. 


(6.58) 


In the second one we use |0^.| < 1 to estimate 

iim Vl [ e'^-F^dX < ||^|| i oclim V(||FPJ| i oo|5 €P (X,)| + |FP|(5 € p(X fc ))) 


(6.59) 


By Lemma 5.4, |F £ |(.B e p(Xfc)) < C 2 §e 2p , and recalling that ||F £a || Loo is bounded, N e < C/e and 
p > 1/2 we conclude 


hmVl [ 0' k {Fl a -F*)ct>dX < ||^limC 26 iV e e 2 P = 0. 
e->0 fc 


(6.60) 


As a result F 9 a F p in and therefore in L°° as e tends to zero. Moreover, by the definition 
of <^e Q , Curl F? a = 0 in \ U/c£> ce (X&). Furthermore, |Curl F^K^X,) X N e N s ne < oo since each of 
the 7V e dislocation points has at most 7V S radial jump sets and each of them has a Burgers vector 
bounded by ne, see Fig. El We recall that the determinant of a 2 by 2 tensor can be expressed 
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as the scalar product of two vectors whose Curl and Div (divergence) are controlled by the Curl of 
the tensor. Indeed, 


det F = (Fn, F 12 ) x (F 2 i, F22) = CF 11 , ^ 12 ) • (^ 22 , —^ 21 )? 


(6.61) 


with Curl (Fh,Fi 2 ) = ^ 12,1 — ^ 11 . 2 ? and Div (F 22 , ~F 2 i) = ^ 22 , 1 -^ 21,2 = Curl (F 2i ,F 22 ). Therefore 
the div-curl lemma (Evans, 1990) implies detFj? a detF p in L°°. 

Next we show that detFj? a —> 1 in L 1 . Since detFl? a = 1 on the union of the cj/, and it is zero 
elsewhere, 

|| detFP, - l|| L i (n /) < \B e p (Xfc)| < TrN e e 2p 0 (6.62) 

k 

since p > 1/2. Therefore detF p = 1 a.e., which is the sought-after result. □ 

In closing we remark that the determinant of the gradient of a SBV function such as (p e can 
be given a distributional interpretation. This formulation is not used in the proof above, but 
illustrates the idea behind the statement detF p = 1. Indeed, in Appendix A we show that if 
< p e G L°° fl SBV(fL;M?) and cj) G C%°(Q), the determinant can be interpreted as 


(det D(p e , (j)) = ~\ [ cof \ 7 ip e D(j)dX + l - f x M • (N x D<j>) dU 1 . 
J O J J 


(6.63) 


((p € ® D(f>) : dcof Dcp e , 


where the precise definition of cp e on J is irrelevant since ipf x [<^ e ] = cp e x [<^ e ]. One can therefore 
choose for <p e (X) on X G J any convex combination of <^+(X) and <^“(X). 

7 Conclusions 


The kinematic assumption F = F e F p pertaining to the modeling of elastoplasticity in the setting 
of large deformations was introduced in the 1960’s based on heuristic arguments and has become 
standard in the continuum mechanics community. However, the lack of a micromechanical justifi¬ 
cation has raised some skcepticism over the years about its validity, in particular questioning the 
existence and uniqueness of the multiplicative decomposition. 

In this paper, we provide for the first time, as far as the authors’ knowledge, a rigorous proof 
of the expression F = F e F p for a general elastoplastic deformation of a single crystal. The proof is 
based on the coarse graining of a mesoscopic description of the deformation, where the dislocations 
and slip surfaces are individually resolved and the displacement field can be treated as continuous 
in all the domain except at the surfaces over which dislocations have glided. At such scale, as 


previously shown by two of the authors (Reina and Conti, 2014), there exists physically-based 
definitions of the different tensors F e , F^ and Fg that are uniquely defined from the microscopic 
deformation gradient cp e and that do not make use of any multiplicative relationship between them 
for their definition. Instead, cp e is described in the compatible regions away from the dislocations 
as the composition of an elastic and plastic deformation, i.e. cp e = o The corresponding 
continuum quantities at the macroscopic scale F, F e and F p are defined respectively as the limit 
of (p e , F e , F^ and F|? as the lattice parameter e tends to zero in the appropriate topology. We show 
that the limiting deformation mapping cp is indeed continuous, that F = F e F p holds in the limit 


30 








with det F p = 1, and that Curl F p represents the dislocation density tensor when it is expressed in 
the reference configuration. 

The results presented in this paper are of great generality, although limited so far to two dimen¬ 
sional deformations due to the mathematical complexity involved in the proofs. In particular, the 
assumptions made consist of standard growth conditions for the elastic energy, dislocations sepa¬ 
rated from each other and from the boundary by a distance of the order of several atomic spacings, 
a finite variation of the elastic deformation and a precise sequence of slips assumed in the compat¬ 
ible regions of the domain. We remark that this latter assumption is only needed for the proof of 
det F p = 1 and it has only been introduced for mathematical simplicity. It is yet to be studied 
in further detail whether or not some of these assumptions are necessary for the final results to hold. 


Appendix A 


Distributional Curl of a second order tensor. Let </> G and F p smooth. Then using 

indicial notation and integration by parts, 


f <t> (Curl F p ) i dX — f (f) e Zkj (F% k dX = - f e 3kj </> >k (F p ) tj dX = - f (D<f> x F p ) i dX, 
Jn Jq Jn Jq 

(•!) 

where is the Levi-Civita symbol. The equality of the first and last term is taken as the definition 
for the distributional Curl. 

Determinant of D(£, with <p G SBV. Let (j) G C£°(fi), and <£> smooth. Then, the weak and 
strong form of the determinant coincide and read 


[ (det D<p)cf>dX = -- [ ^ T (cof D<p)D(f>dX = -- [ (ip <g> D<j>) : (cof Dtp) dX. (.2) 
Jn 2 Jq 2 Jq 

Indeed, 

P (cof D V ) D</> = ( ifii ) ( Fplz ) ( % ) = ( ^ V* ) 

= PlP2,2(f>,l — PlP2,l4>,2 — P>2 P 1,2^,1 + P2Pl,l(f ) ,2- 

(.3) 


^2,20,1 — P2,1<P,2 
— ^1,20,1 + Pi, 10,2 


Therefore, after integrating by parts, 

- [ p T (coi Dip)D<f>dX 
Jn 

= / (Pl,lP2,2 + PlP2,21 ~ <Pl,2<P2,l ~ PlP2,12 ~ P2,lPl,2 ~ P2Pl,21 + P2,2Pl,l + P2Pl,n) <P dX 
Jn 

= 2 / (pi } iP2,2 - Pi,2P2,i) <fidX = 2 / (det D<p)<j>dX. 

Jn Jn 


(•4) 
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For p e SBV, cof Dip = cof \7pC 2 + cof(fo>] <g> N) U l [j, with 


~Mm 

IflijNi 


Mm Mm 

lV2 

-M -m 


cot ([^?] <g) n ) = cof ( ri;; 1 y 2 
' 1 Mm mn 2 


[<P1 


m 


Mm 

-I^i]iV 2 


= l‘pi ± ^N- 1 , 




cof (Ipj ®N)D4>=(p- M X ) (n X • D<j) 


vi W -M ^ (( -m 
V2 ) ' V M )) Vv m 

= (~V X M) • (N X D4 >), 


therefore giving (6.63). 



(.5) 


(. 6 ) 
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